m 
o 

C 



o 

c3 



i 

c 

o 
o 



> 

(N 

oo 
m 
p 

o 



x 

S3 



Photo-mechanical energy conversion using polymer brush dissociation 

J. M. Deutsch 1 

Department of Physics, University of California, Santa Cruz CA 95064 

A device is investigated that continuously and directly converts light into mechanical energy, using polymers and pho- 
todissociation. A polymer brush tethered to a surface, is brought into contact with a parallel plate a small distance above 
it that contains reaction sites where photodissociation of bound polymer and light can occur. Under the appropriate con - 
ditions, the collective effect of these polymers is to apply a force parallel to the plates, converting incoming light into 
mechanical work. Numerical work is carried out to understand this effect, a three dimensional Langevin simulation, 
solution to the Fokker Planck equation, and a one dimensional Monte Carlo simulation. Theoretical analysis of the 
Fokker Planck equation is used to study a model where equilibration of the unbound state occurs and equilibration to a 
metastable equilibrium is achieved in the bound state. It is shown that the work per cycle can be made much larger than 
the thermal energy but at the expense of requiring a greatly diminished photodissociation rate. Parameters are discussed 
in order optimize mechanical energy conversion. 



I. INTRODUCTION 

There are many proposals to convert light into mechanical 
energy using smart polymeric photo-responsive materials 1 " 5 
or the synthesis of individual molecules that lead to rotary or 
linear motion 6 . 

Photoresponsive materials use a conformational change of 
a macroscopic collection of polymers in response to light or 
an external change in temperature and pH. This will cause a 
change in volume. If this process is reversible, which it is in 
many instances, this device can be cycled to derive mechanical 
energy^ - —. 

In general, there are many scenarios leading to a change 
in configuration of a molecule in response to light 2 : cis- 
trans isomerization, zwitter ion formation, radical formation, 
ionic dissociation, and ring formation/cleavage. If such light- 
sensitive elements are incorporated into a macromolecular 
system, this can in principle, through thermodynamic cycles, 
to conversion of light to mechanical power. 

Remarkable advances in the design of molecules have 
led to prototypes for light-driven synthetic molecular mo- 
tors^. Many of these are based on photoisomerization reac- 
tions that enable rotary motion of molecules when combined 
with thermal rotation steps 6-9 . It has also been possible to 
make threading-dethreading systems 610 , and shuttle molecu- 
lar rings&i!. There has even been progress on a macroscopic 
level, such as creating droplet motion by modifying surface 
properties 12 or inducing mechanical deformation 13 similar to 
what is achieved with photo-responsive materials^ - — . 

Much of the research in this area is inspired by biologi- 
cal machines such as myosin II that work by quite a different 
principle than man-made macroscopic motors 14 . Biological 
motors use chemical energy rather than photons to produce 
mechanical power however this difference does not affect the 
basic mechanism of operation. A myosin head can be thought 
of as being in two states, bound or unbound. Thermal noise in 
the unbound state can cause the head to bind to actin, produc- 
ing a force. The hydrolysis of ATP releases energy causing 
the head to return to the unbound state. The biochemistry of 
a real motor protein is considerably more complicated, but by 
simplifying this description to one involving only these two 
states, the motion can be analyzed, and it is easily seen that 



electromagnetic energy can be used instead of chemical en- 
ergy^. 

In this paper, I investigate the use of photons in powering 
a two state motor system similar to biological motors. I con- 
sider creating motion between two surfaces that are very close 
together by placing an asymmetric polymer brush between 
them. This could potentially have advantages. A continu- 
ous source of light would create a motion that is constant on 
a macroscopic scale, that could for example, rotate the sur- 
faces relative to each other. It would not require any addi- 
tional mechanisms to keep it moving, other than the micro- 
scopic motion of molecules between the plates. 

The device proposed here falls into a distinct category dif- 
ferent than the experimental approaches above. It is not 
a macroscopic smart material that changes properties in re- 
sponse to external stimuli. The device described in the next 
section works on the scale of an individual polymers, and the 
force generated on the surfaces is the sum of these molecules 
acting independently of each other. However it is unlike the 
chemical synthesis approach that requires different states of 
isomerization. The principles that it relies on are robust and 
just require photodissociation of polymer to a binding site, 
and, as in the case of biological motors, some degree of asym- 
metry. This asymmetry, plus a disruption of thermal equilib- 
rium due to energy input are the two factors needed to convert 
energy from chemical or electromagnetic energy, to mechani- 
cal 15 . 

The paper is organized as follows. In Sec. [II] the photo- 
mechanical system is described and rough estimates of its op- 
eration are given in Sec. Hm including a discussion of its po- 
tential efficiency. In order to illustrate the characteristics that 
need to be understood, a three dimensional simulation of this 
device is carried out In Sec. [TV] At this point the physics of 
this system is investigated in more detail, starting with general 
considerations, in Sec. |V]of its steady state behavior using 
Fokker Planck equations. The results will be used in subse- 
quent sections. An exact description of this system is possible 
in steady state and the resulting differential equation is solved 
numerically in one dimension in Sec. [VI] This allows to un- 
derstand better how effective the asymmetry in the force pro- 
duces power, and influences the direction that is taken in the 
rest of this work. Sec. IVIII studies a useful limit that is hard 
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to probe numerically, but allows to understand how the effi- 
ciency of this system is related to spring stiffness, metastabil- 
ity, relaxation times, and the photo-dissociation rate. We then 
use this model to understand the efficiency, Sec. IVII Al using 
a one dimensional Monte Carlo model that is in this regime. 
Finally in Sec. I Villi we conclude on how the physics that has 
been learned might be useful in optimizing experimental pa- 
rameters for such a system. 

II. THE SYSTEM 

The components needed to do constant photomechanical 
energy conversion are are illustrated in Fig. [T] 

• A flat plate of semi-flexible polymer brushes. Poly- 
mer brushes are polymers, each with one end tethered 
to a wall in a suitable solvent. 

• A parallel plate right above the brush. These poly- 
mers are put in contact with a parallel plate close to the 
surface so that their ends are able to interact with it. 

• The parallel plate contains an array of photoreac- 
tive binding sites. It is crucial that these polymers bind 
with the surface in an asymmetric way, so that the av- 
erage binding orientation of each is the same, and not 
perpendicular to the plate. 

• At least one of the plates is transparent. Light causes 
the unbinding of polymer ends from the photoreactive 
binding sites. 

• Binding catalyst. To control the rate at which binding 
occurs, the binding of the end of the polymer to a bind- 
ing site can be facilitated by the use of a catalyst. The 
concentration of the catalyst can be used to control the 
rate of binding. 

In most of the discussion below, the polymers can be con- 
sidered to be separated from each other by a sufficient distance 
so that we can ignore inter-chain interactions. The absorption 
of light causing unbinding will happen asynchronously. The 
ends tethered to the lower plates are at positions that are either 
random or incommensurate with the binding sites of the upper 
plate. Altogether this means that the motion of each polymer 
is uncorrelated with the others in the system. 

The total effect of the forces acting on the plates can be used 
to perform work against a force acting on the upper plate. Be- 
cause we assume that the number of polymers contributing is 
very large and are asynchronous, the net velocity of the upper 
plate v will be constant, as will be the net force acting on the 
plate. 

In the following we will analyze the system and model it in 
different ways to better understand the power generation. 

III. ESTIMATE OF SYSTEM PARAMETERS 

We will now make a crude estimate with realistic parame- 
ters, how well this device should work. We envisage separate 



one dimensional tracks of the kind shown in Fig. Q2b). The 
spacing between the tracks could be greater than the distances 
between polymers on a single track. In reality, it might prove 
more efficacious to manufacture a square array of polymers 
and binding sites but we will allow two separate length scales 
in our analysis below. 

We will assume a solar intensity of / = 600W /m 2 , and an 
average photon energy of e = 2eV. We will take the relax- 
ation time of a polymer to be T = 10 s, which corresponds 
roughly to a polymer size of 3nm. The number of unbind- 
ing events, if every photon was absorbed at a binding site, is 
//e«2x 10 21 / (m 2 s). If the density of polymers is a, then 
for all of these events to be utilized requires I/e — ct/t, or 
a = 2 x 10 14 /m 2 , which is a separation of \j \fo sa lOnm. 
The amount of energy per step that is gained by a photo- 
dissociation event is of order IcbT. With such parameters, this 
is expected to yield an efficiency of order 1%. However, as 
we will show below, longer relaxation times allow for more 
energy per step, and this requires a higher polymer density. 
However the two dimensional density is limited by inter-chain 
interactions. In the case considered here, this density could be 
increased roughly by three orders of magnitude. However the 
relaxation time for this system depends exponentially on the 
force generated, so a three order of magnitude increase in re- 
laxation time would only increase the efficiency by a factor 
of about 7 furthermore, the detailed analysis presented below 
suggests that with optimal conditions, the efficiency of con- 
version is about 7%. To increase the efficiency of solar con- 
version further might require stacking devices. 

There are other means of increasing the solar conversion 
efficiency by first converting solar radiation to lower energy 
photons. For example, two methods for doing this are fluo- 
rescent down-conversion of solar photons 16 , or thermophoto- 
voltaic down-conversioraii 

Further speculation on device efficiency is premature as 
there will undoubtedly be many unforeseen technical prob- 
lems that will likely provide other obstacles to increasing 
device efficiency. However an efficiency for direct photo- 
mechanical conversion of about 7% is still a useful amount 
of power comparable to photovoltaic conversion with amor- 
phous solar cells, and also because power is lost in electro- 
mechanical conversion, which is not a problem with direct 
mechanical energy conversion. 



IV. THREE DIMENSIONAL MODEL 

We start by simulating a three dimensional model of this 
system. There are two components to the system, the polymer 
and the binding sites. The polymer chain is modeled as having 
N links of fixed length, and is semiflexible with chain stiffness 
K. Denote the coordinates of the ith bead as r,. The elastic 
potential for the middle of the chain for the ith bead is 

U E (i) = -|(|r i+2 -r ; | 2 + |r,_ 2 -r,| 2 ) (1) 
where C is the elastic constant. 
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FIG. 1. (a) A schematic of the light energy to mechanical converter. A polymer brush of semiflexible chains is put in contact with a plate a 
height h above the base, containing an array of anisotropic binding sites. Polymer ends are trapped by unlikely thermal fluctuations, applying 
a non-zero average force parallel to the plate. Photodissociation of the polymer ends with the binding sites releases the chain end that will 
move until it finds another site. The upper plate moves at a velocity v, thereby generating power. A single track is shown, and the separation 
between binding sites is L. (b) A three dimensional sketch of the binding sites. The binding of the polymer ends (not shown) takes place on 
the oval surfaces protruding from the plate, and do so in an asymmetric manner. The sites are organized into one dimensional tracks, where 
the distance between tracks need not be the same as L in (a). 



The binding potential of the polymer has two components, 
an isotropic component £/, and directional component Ud- 
Ui(r) is short range with a length scale r s and scale V a 



Ui(r) = -¥ztf-r 



„2\4 



(2) 



and Ud uses a direction n a so that the difference between the 
last two end beads Ar = ri ro will give a minimum in Ud 
along that direction 



U d (r ,ri) = ( 



Ar n fl 
I Ar - n„ 



(3) 



The end attached to the lower surface is always bound. The 
other end can bind to a periodic linear array of binding sites 
equally spaced at a distance L. 

There are two states that the system can be in, unbound, 
0, and bound 1 . There are two parameters that control bind- 
ing, the rate at which dissociation occurs co, and the rate that 
dissociated ends can be rebound c\. 

The distance between the lower and upper surface is h and 
their relative velocity is v. 

The model was simulated with a Langevin equation, at fi- 
nite temperature T. Although inertial effects are typically 
small at these microscopic scales, it was included for com- 
pleteness. An algorithm was used that efficiently updates this 
system with fixed link lengths^. 

The simulation was run to determine how the force and 
power generated are influenced by the plate velocity. Fig. [2] 
shows the results of four simulations with different values of 
Co and c\. The parameters used were: 4 links, h = 1, spac- 
ing between binding sites L = 2, the direction of n„, is 7t/4, 
V a = 30, link length of 1, particle mass of 1, C — 1 .5, r s = 2.0, 
and a coefficient of damping of 10. 

The force generated is highest approaching zero velocity 
and decreases until it becomes negative, at which point it 



is taking mechanical energy to move it at such high speeds. 
Qualitatively, this is the point where frictional drag dominates 
over photo-energy conversion. 

Note that the highest power P = /v is seen for co = c\ = 
0.05. The force extrapolated to v = is / = 2.07 and is about 
three times less than the force seen in Fig. |2jb) where cq — 
0.05 and co =0.01. 

Higher maxima in the power are generally seen with in- 
creasing co and c\ . As an example, the simulation was run for 
small elastic constant C = 0.5, at different values of co and c\, 
but with all other parameters the same as above. This higher 
power is still at the expense of efficiency, as higher dissocia- 
tion rates imply a higher photon flux. 



V. STEADY STATE PROBABILITY DISTRIBUTION 

To understand the above simulation in more detail, requires 
a better understanding of the mechanisms involved. There 
have been a very large number of works on the theory of mo- 
tor proteins 1419 , and we will follow the two state approach 
mentioned above of Prost et. a/ . 15 ' 20 . The major difference is 
that they considered the two states to be periodic potentials, 
whereas here we model the system to more closely mimic the 
particular device we are investigating allowing the exploration 
of force versus velocity. Hence instead of considering the po- 
sition of a particle moving between two periodic potentials, 
we consider the unbound state to have a free end described 
as moving in a bound potential around the tether point. Like- 
wise, because of this tethering, the bound state potential is not 
periodic but has a periodic component as we will describe in 
detail below. Here we consider the probability distribution for 
this system in steady state, which is described by a Fokker 
Planck equation. 

Initially to calculate the power that is produced, we will 
concentrate on the low velocity limit, with the two plates mov- 
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Power 



0.0000 0.0005 0.0010 0.0015 0.0020 0.0025 
V 



(b) 



Force 
Power 



01)00 0.002 0.004 0.006 0.008 0.010 0.012 
V 




(c) 



000 0.002 0.004 0.006 0.008 0.010 
V 



(d) 




0.0000 0.0002 0.0004 0.0006 0.0008 0.0010 0.0012 
V 



FIG. 2. The average force / and the power P = /v measured as a function of the relative velocity between the plates v. For clarity, the power 
is multiplied by 2000. There were 4 links, with a plate separation of h = 1, spacing between binding sites L = 2, the direction of n„, that is, 
the binding angle is 7t/4, V a = 30, link length of 1, and elastic constant C = 1.5, potential range r s = 2.0, and a coefficient of damping of 10. 
(a) the rate of dissociation c = 0.01, rate of rebinding cj = 0.05. (b) c = 0.05, c\ = 0.01 (c) c = 0.05, ci = 0.05 (d) c = 0.005, ci = 0.05 




(a) 



0X101 0.002 0.003 0.004 0.005 0.006 
V 




(b) 



000 0.005 0.010 0.015 0.020 0.025 0.030 0.035 
V 



FIG. 3. The average force / and the power P = fv measured as a function of the relative velocity between the plates v. The same parameters 
and scaling are used as in Fig. |2]except C = 0.5 and (a) co = c\ = 0.05, and (b) co = c\ = 0.4. Higher maximum power is seen for larger 
unbinding rates. 



ing so slowly that this motion does not affect chain conforma- 
tions appreciably. In this way, we can consider the average 
force exerted between the plates by a polymer in steady state 
in the limit v = so that the tether point is not moving. It is 
most convenient to let the point at which the polymer is teth- 



ered to the bottom plate be a variable parameter r' = xfx. With 
this point fixed, we consider the distribution of the other end 
of the polymer r. We will assume that the internal dynamics 
of the chain chain are much faster than the binding and un- 
binding rates, so that the only degrees of freedom are r, and 
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if the polymer is bound, ; = is unbound and i = 1 is bound. 
Therefore the probability distribution of the system can be de- 
scribed by a function P s (r;r'). The equations describing this 



hence 



an 



.15 



d t P (r) = V • (V - fo)Po(r) - c P (r) + Cl P x (r) (4a) 
d,P x (r) = V • (V - f i )P 1 (r) + c P (r) - cifl (r) (4b) 

where the last r' argument of P has been left out for notational 
simplicity. The units here absorb the diffusion coefficient D 
together with the time, that is t as used here and below is really 
D times the time. The force terms have absorbed a tempera- 
ture factor T, that is f,, is really the force times 1 /ksT . Those 
forces, fo and f i are the total forces acting on the upper end 
of the chain. When the system it is unbound is the force of a 
(possibly) nonlinear spring 



fo=f,(r-r'; 



(5) 



When the system is bound we have total force is the sum of 
the spring force and a periodic force representing the binding 
potential 



f 1= f,(r-r')+f p (r) 



(6) 



These forces are assumed to be conservative, f s and f p are 
derived from potentials that respectively are V s and V p . 

In steady state the left hand sides of these equations are 
zero. We can eliminate the last two terms on the right hand 
side to obtain. 



where 



O0P0 + O1P1 =0 



d/ = V-(V-fi), z = 0,l. 



(7) 



(8) 



We can now calculate the average force exerted by the two 
potentials, is zero. In steady state, we multiply Eq. |7] by x 
and integrate with respect to x, y and z. Then using integration 
by parts, we have boundary terms at infinity. Because we are 
assuming that the spring potential grows without bounds, this 
confines Pq, and Pi, to a neighborhood around r', so that the 
boundary terms vanish. We are then left with 



f>= /(f Po + fiPiK / r = 



(9) 



as is expected because r is confined to a region of space so that 
the average velocity, and hence average force, will be zero. 

We can also determine the total fraction of time spent in 
the bound or unbound states in steady state. First, because the 
probability of being in any state is unity, 



(P (r)+Pi(r)K / r=l. 



(10) 



Then by integrating Eq. |4a]over all space, the derivative term 
integrates to 0, giving 

j(-c Q P Q (r)+c l P 1 (r))d d r = (11) 



P (r)d d r = — — , f P x {r)d d r ■■ 

CQ+Cl J 



CO 



CQ + Cl 



(12) 



To obtain the power produced by this device, we are not 
interested in the total force acting on the upper chain end be- 
cause this includes the binding potential, but the average force 
due to the spring acting on the lower plate (f s ). We would like 
to calculate the work done in moving the lower point r' by one 
period of f p . After moving one period the system is statisti- 
cally identical to its starting point, and this method can there- 
fore give the work done in moving n such periods. Denoting 
the period of i p by L, we would like to calculate 

W,= t t-(t s )dx! 



x.-f s (r,x'x)(P (r,x'x)+P l (r,x'x))d 3 rdx' (13) 



and using Eq. [9] 



Wi = -J J x f p (r)Pi (r,yx)(f J rJi'. (14) 

In thermal equilibrium, where the transition rates cq and 
Ci are both zero, we recover the Gibbs distribution. Let us 
assume, that the system starts, and therefore remains in state 
i = 1 . Then 



P (r,r')d 3 r=— d\ 
where the partition function 

Z(r')= / e- Vl d 3 r 



(15) 



(16) 



which will also have periodicity of L. In this case Wl can 
be easily calculated because Pq = and f s = — W s (r — r') = 
V'V s .(r-r / ). So 



W L = jJ(d,V s (r, X %) W) 

= -r(lo g Z(L)-logZ(0))=0 



d\dx' 



(17) 



as it must be by the second law of thermodynamics. 



VI. ONE DIMENSIONAL SOLUTION 

In order to investigate how the power conversion depends 
on the forces acting on this system, it is important to simplify 
the three dimensional model to obtain a minimal model that 
depends on far fewer parameters. Therefore we investigate 
this model in one dimension. 

A key point to understand is how asymmetry in the form of 
the force produces power. With symmetric spring and binding 
potentials, it is easily seem by symmetry, that no net power 
can be produced from this system. We now ask how asym- 
metry affects the results. We will see that even with a large 
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asymmetry in the spring potential, the work defined by Eq. 
n~3lis very small. Simulations using Monte Carlo or Langevin 
equations are too noisy to provide good estimates. We there- 
fore instead use a more analytical approach. 

The coupled Focker Planck Eqs. |4](a) and (b) can be solved 
to produce an equation only involving one distribution func- 
tion, Pq in steady state. Using the definitions in Eq. [8] we can 
eliminate Pi . 



(d 1 do-c 1 d -cod 1 )p = o 



(18) 



which is a fourth order equation in spatial variables. 

Now we restrict the analysis to one dimension. In this case, 
Oi = d x (d x — fi) so we can integrate with respect to x. We note 
that because the spring confines Pq to a localized region, it will 
go to zero as x —> ±°°. Therefore the integration constant must 
also be zero 

((<?,- - fi)d x (d x - /o) -a(d x - /o) - co^-Zi))^ -0(19) 

which is a third order linear differential equation. 

Eq. [I9]was solved by the shooting method 2 1 . The boundary 
conditions were obtained by considering the system far from 
x' where Pq is very small. In that domain, f p was artificially 
cutoff so that /o = f\ = f s . Because the potential there is no 
longer changing between the two states, the solution is that of 
a system in thermal equilibrium. The solutions were required 
to match to these thermal solutions in this regime, far from 
x'. The equation was solved with three different initial condi- 
tions. An appropriate linear combination of these were con- 
structed to match the boundary conditions as just described. 

The periodic binding potential that was used is 



V P =A\ cos(x) — A2sin(2x) 



(20) 



We first consider asymmetry in V s but with symmetric func- 
tions V p , that is, A 2 = in Eq. [20] 




FIG. 4. An asymmetric spring potential used to tether the chain. It 
has a second dip at approximately x = 2 



Q- 0.6 




Vs{x) = \ kx2 -i+{ b {x-d)r 



(21) 



The first term describes a linear spring with spring constant 
k, the second adds an asymmetric dip. The parameters are 
chosen so that this dip is close in potential to the one created 
predominantly by the linear term, k = 4, a = 8, b — 2, and 
d = 2. The function is plotted in Fig. H 

Plots are shown in Fig. of P\ (x) for four values of x' 
within a period. 

By integrating using these distributions, Eq. [14] the work 
can be obtained. With cq = 0.025 and c\ = 0.05 the work 
W L = 0.00080954. With the c/'s 10 times those values, c = 
0.25, and a = 0.5, W L = 0.0025. 

Now we consider the case of a linear spring so that the non- 
linear parameter a = in Eq. |2T] Instead we make the periodic 
potential asymmetric by setting A\ = A 2 = 2 in Eq. [20] still 
with co = 0.25, and ci = 0.5. Fig. [6]plots Pi (x) for four values 
of x 1 within a period. In this case, the work Wl = —0.1596. 

What the numerical results have shown is that asymmetry 
in the spring potential is quite ineffective at producing work, 



FIG. 5. Plots of the probability distribution of P\ as a function of 
position x for different values of x' for the asymmetric spring model 
Eq.M 



whereas asymmetry in the periodic binding potential is much 
more effective. We shall use this result below and concentrate 
on systems with symmetric spring potentials, but asymmetric 
potentials. 



VII. UNBOUND EQUILIBRATION MODEL 

It is worthwhile to understand in more detail, what con- 
strains the maximum force of photo-mechanical conversion 
by tuning the potentials employed and the binding rate. And 
to ask how the design will depend on the flux of photons. We 
are limited in our choice of the potentials Vq and V\ that both 
must be bounded. In this system, the unbound potential is not 
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FIG. 6. Plots of the probability distribution of P\ as a function 
of position x for different values of x' for an asymmetric periodic 
potential, but a linear spring. 



periodic which limits the amount of power that can be gen- 
erated. In cases considered earlier 15 , where both potentials 
are periodic, it is possible to get much more efficient motion 
by alternating between states with different potential maxima. 
However this is not relevant to our system. 

To understand this better, it is useful to examine a limit 
where we can treat the system analytically. We therefore ex- 
amining the case where the binding rate of unbound chain is 
sufficiently that we can regard it in thermal equilibrium. 

The model is illustrated in Fig. [7] We take the unbound 
potential to be that of a linear spring below a cutoff 

V s =Vo = P /2 ' tf -<°<*<k (22) 
I °o j otherwise. 

we choose the spring coefficient k and the cutoff length l x such 
that l c exp(-klc /2) < 1. Below, we will take l c = L/2 or oo. 

The periodic potential V p is taken to be localized at periodic 
points (with a separation of L) that rapidly vary from a large 
maximum y mal to a minimum V,,„„, as shown. We take the re- 
gion over which this happens to be negligibly small compared 
to other length scales in the problem. 

To understand how equilibration can occur in this model, 
we can rewrite Eqs. |4]in the steady state limit 

V • (V - fo)Po(r) - c P Q (r) = -aft (r) (23a) 

V • (V - f j )P 1 (r) - ciPi (r) = -c P (r) (23b) 

Although this is time independent, we consider the related 
time dependent equations for variables P, (r,f ) 

d,P (r) - V • (V - f )/b (r) = (24a) 
o »,A(r)-V.(V-f 1 )Pi(r)=0 (24b) 

where P,, i = 0,1, is the Laplace transform of P,, Pi(s) = 
J^{Pj}, where the conjugate variable for equations (a) and (b) 



FIG. 7. Illustration of the kind of potentials employed in the Un- 
bound Equilibration Model. The periodic potential V p is nonzero 
only on small regions on the x-axis. It shows a large peak of 
height Vmax an d large negative dip V m ,„ . There spring potential Vb is 
parabolic, except that it has a cutoff where it becomes infinite when 
stretched by more than Z c . The corresponding probability distribution 
Pq (grey curve) is assumed to have relaxed to equilibrium. The spring 
potential is shown at two different position, when it is centered at x' 
and L — x'. 



are s — cq and s — c\, respectively. The initial conditions on 
each equation are 

P (r,t = 0)=c 1 P 1 (r,s = c 1 ) (25a) 
P l {r,t = 0)=c Q P Q {r,s = c ) (25b) 

These equation have a direct physical interpretation. The 
left hand sides in Eqs. |24]are those of particles in diffusing in 
potentials but with conservation of particles. For long enough 
times, independent of initial conditions, the solution to these 
equations will go to thermal equilibrium given by the Gibbs 
distribution. Eq. I24al describes diffusion in a quadratic poten- 
tial. We require that we are probing this at long enough times 
t, so that it will have nearly reached this equilibrium state. 
Call this longest relaxation time To- The solution to the diffu- 
sion equation P(r,t) for either of Eqs. [24] can be written as a 
sum over spatial eigenfunctions, <j>„ (r) that decay at different 
rates A„, (arranged to be monotonically increasing): 

P{r,t) = f>„(r)exp(-A„r), (26) 

of which the smallest A, Ao = 0, corresponds to the equilib- 
rium state. The relaxation time is Tq = Thus the Laplace 
transformed variable P can be written 

'M=E^)^7-. (27) 

If the n — term is to dominate, we therefore require s = 
c o "C 1/ lb- The physical interpretation of this condition is 
that binding typically occurs only after many relaxation times 
of the unbound end. 

Pq(t,s = co) will be dominated by the first term in Eq. 127] 
and therefore proportional to the eigenfunction 0o( r ) which 
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is the equilibrium distribution and is exp(— Vb). An impor- 
tant simplifying point in the above approach is that the precise 
form of the initial condition Eq. I25al is not not important, but 
because of particle conservation, only the total area under P\ 
affects the result. 

Now we consider the solution for Pi . If we consider un- 
binding times much longer than the relaxation time in the 
bound state, we also arrive a thermal equilibrium, which as 
was shown by Eq. [T7] to lead to no work being performed. 
Instead we will consider situations where there are very long 
lived metastable states. In Fig. [7] the potential seen by a par- 
ticle in Eq. [24b] is V\ = Vo + V p . If a particle starts between 
and L, it will remain trapped in that region for a Kramer's 
time which (ignoring algebraic pre-factors) depends on xf, but 
has a minimum value of %„ oc exp(V„ !M ). By choosing large 
enough V max this can be made arbitrarily long. 

A particle in such a metastable state will relax to a 
metastable equilibrium, obeying Eq. [26] that will eventually 
fail for times t > X m , In this expansion, the longest relaxation 
time to this metastable state Ti = 1 /X\ , will be taken to be 
much smaller than x m . Thus the above argument on the range 
of co can be used mutatis mutandis to restrict the unbinding 
rate to l/r m <ci < 1/ti. 

In this regime we can understand the solution to Eqs. 124b I 
and I25bl by considering the corresponding Greens function 
G(r,f;ro). We replace the initial condition Eq. |25b| by 

Gi(r,* = 0;r ) = S(r-r ). (28) 

We can then obtain P\ from the Greens function from 

A(r,r) = j Gi(r,t;r )c Po(ro,s = c )d d r Q . (29) 

For a wide range of times, and regions of r', the G\ will 
approach the same metastable state. A particle starting at any 
point within a certain region will end up stuck in the same 
state and hence approach the same metastable equilibrium. 
Confining our attention to the one dimension model of Fig. 
[7] for < xq < L/2 the solution will be strongly localized at 
the minimum x = 0. For L/2 < xq < L the effects of V„ are 
negligible and P\(x) exp(— Vo). This is because the poten- 
tial Vo in Eq. [22] is cutoff which does not allow the particle 
to visit the x — region. Hence the effect of V p is only seen 
close to x = L where it provides a strong repulsion, but over 
a negligibly small region of x. For —L/2 < xq < 0, V p also 
does not contribute. For L < xq < 3L/2 the solution will be 
strongly localized at the minimum at x = L. Because of parti- 
cle conservation, the area under G\ is always unity. 

A physical interpretation of the above equations in terms 
of a one dimensional one particle system can now be made 
using the Laplace transformed variables and the metastable 
limit considered above is now apparent. In the unbound state, 
the particle reaches thermal equilibrium relaxing to the Gibbs 
distribution, Pq °= exp(— Vo). Then the periodic potential V p 
is suddenly added in. The position at the time of switching 
is labeled xq. Depending on which interval xq is in, x will 
equilibrate to the corresponding metastable equilibrium with 
P\ « exp(— Vi ) and is completely confined to that interval. The 



relative probabilities being in one of the three above regions 
is obtained by the area under Pq for that interval. 

Now we know how to determine Pq and Pi , we would like 
to calculate the work Wl given by Eq. Q~3] Because of the 
symmetric form assumed for f s , the /qPq term in the integrand 
gives zero contribution and we are left with 



W ' = Jo J M x ~ x ') p dx,x')dxd£ (30) 
and consider first how to calculate the inner integral 



f,(x') = J / (*-jOPi(VM* 



(31) 



where x 1 is the position of the tethered end. Eq. [22] gives 
fo( x ) — —kx for |x| < L/2. 

Using the prescription we have found for Pi (x^x 1 ), which is 
simplified by the above physical interpretation, we can parti- 
tion this integral into the different x-intervals of metastability: 
/_ = [-L/2,0], Iq = [0,L], and /+ = [L.3L/2]. The value of 
P\(x,x') depends on the probability of initially being trapped 
in one of those three intervals. Because of the cutoff we have 
imposed on Vo, only two intervals need be considered for a 
given value of x 1 . For < x' < L/2, only intervals /_ and Iq 
occur. The probability that i£L given x' is 

E (xf) = Prob(x € /- \x) — j po(x — x r )dx — J po(x)dx 

(32) 

and the probability that x e Iq that is, P(x E I \x') = \-E(x'). 
Here po(x) is proportional to Po(x) but normalized to unity, 
to simplify the presentation. We can obtain the values for 
L/2 < x 1 < L by symmetry, so that P(x G I+\x!) = E{L — x!) 
and P(x £ Iq\x') = 1 — E(L — xf), E(x) is simply related to 
the complementary error function in the limit considered here, 
where the effects of the cutoff in the potential will have negli- 
gible effect, E(x) = ^erfc{^/kx). 

There is a symmetry in many of the quantities considered, 
as shown in Fig. [7]where the potential Vo and corresponding 
probability distribution po is for the tethering point at xf and 
at L — xf. Therefore it is convenient to consider // (xf) +// (L — 
x 1 ). This can be written as 



kCQ 

co + ci 



■(-E(x'){x)'+E(x')x' 
-(l-E(x'))x' + (l-E(x'))(x)") (33) 



where 



f_ a ,(x-x')p (x-x')dx = j"xp {x)dx (34) 



J-ooPo(x — x')dx 



E(x') 



and 



£jx + xf)p (x+xf)dx -E(x') 

(x) = t: = 7—rrW ■ (35) 

j\ P Q(x + x')dx 1 -£(*') W 

The four terms in Eq. |33lcorrespond to contributions from 
respectively regions /_, I + , Iq, and Iq. The first and third term 
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are contributions from fi(x') and the others are from fi(L — 
x'). The factors involving c give the correct normalization 
according to Eq. [12] Combining the above equations, 



2kc f L ' 2 



CQ + Cl 



/ 2E (x) ((x) — x) + x'dx . 
Jo 



(36) 



In the limit of large L, which we are considering by virtue 
of the condition on the cutoff in Vq imposed by Eq. [22] the 
integrand simplifies because (x)' becomes exponentially close 
to x 1 , and the only term remaining is x 1 . Thus for large L, 
Wl = (co/(cq + c\ ))kL 2 /4 The factors involving the c's repre- 
sent the fraction of time spent in the bound configuration. The 
last term increases quadratically with L. This result is mis- 
leading if not taken with the appropriate limits that have been 
assumed in its derivation. The factor (cq/(cq + ci)) is very 
close to unity as we are assuming that the relaxation time in 
the unbound state is much faster than in the bound state, hence 
co c\ . However the work Wl to move a distance L was as- 
sumed to be in the adiabatic limit, and here the time scales 
associated with bound state relaxation are exponentially long. 
This is because to reach this metastable equilibrium the parti- 
cle has to hop over barriers of size Vq(x'), see Fig. [7] Hence 
the longest relaxation time for this is at x' = L/2 and is of 
order exp(fcL 2 /4). Note that this does not contradict our as- 
sumption that we are still in a region of metastability, which 
requires times much less than exp(V max )- But for this to work, 
we require that V max ^>Vq(L/2) 

Now we consider the case where the spring length cutoff 
l c —> °°. The disadvantage of this is that the chain end can, 
in principle, hop over many barriers ending up arbitrarily far 
from the tether point, and that these hopping times should be 
included in the above analysis. However, the probability of 
such a hop becomes negligible when the probability of find- 
ing the chain end there is small. Hence we still have a clear 
separation of time scales between metastable states as dis- 
cussed above, and fully equilibrated system, which requires 
surmounting the energy barrier V majc . In this case, we can 
therefore assume that when the chain end x' is between nL 
and (n + 1)L, it will strongly localized atx' = nL. Therefore 



f l (x') = -k A n (x)(nL~x) 



(37) 



where A„(x') is the probability of initially finding the chain 
end between nL and (n+1 )L, 

A n (x)=E(nL-x)-E((n + l)L-x). (38) 

It is easily seen that A_„(L — x) — A„(x). Using this and Eq. 
[37] 

oo 

f l (L-x) = -k(-f l {x)-L £ A n (x')) = -f,(x)+kL. (39) 

« = — oo 

To obtain the work, we follow the same procedure as above 
and consider // (x 1 ) + fi (L — x 1 ), which here is just kL. There- 
fore in this case, 



W L = 



co kL 
co + ci 2 



which in this model exactly, for all k and L, given that we are 
in a region of strong metastability. 



A. Efficiency in large power stroke limit 

We now are in a position to answer a central question about 
the performance of this kind of device: is the efficiency lim- 
ited by the small value of kgT compared to photon energies? 
We have seen from the above analysis that it is possible to 
get arbitrarily large forces developing at the expense of ex- 
ponentially slow operation. However this is in the limit of 
infinitesimal plate velocity v. In contrast, the power obtained 
is instead average force times this, (/)v and we would like to 
operate the device at the velocity of maximum power, which 
necessitates the operation of it far from equilibrium, because 
compensating the increase in v is a decrease in (/) due to dis- 
sipation. It could be that the optimum velocity of operation 
decreases very quickly with kL 2 /2 meaning that the device 
becomes increasingly inefficient as the spring constant k is in- 
creased. This would make it impossible to harvest more than 
of order ksT energy per cycle. 

The efficiency is defined as the ratio of the amount of power 
produced to the amount of power put in. The amount of en- 
ergy needed to dissociate a polymer end from a binding site 
is V m i n . Binding to V, m „ must produce an energy less than 
that of an unbound polymer. The power put in is V ml „c where 
1/c = 1 /cq + \jc\. Therefore the efficiency is 



ri = 



C Vmin 



(41) 



To investigate this problem further, the one dimensional 
model of the last section with l c — » °° was implemented us- 
ing Metropolis Monte Carlo. We chose the periodic potential 
V p to vary as 



V P {x) 



-4V™„f(l-f), 
-4V mfl 4(l + f), 



0. 



if <x < 8. 
if >x > -8. 
otherwise. 



(42) 



(40) 



Here we set 8 = 0.1, Vmax = Vmin = 100, and L = 1, In order 
to preserve diffusional dynamics, steps in x were attempted 
uniformly in the range [—0.025,0.025] ensuring that the pe- 
riodic potential cannot be jumped across in one move. One 
move increased the time by .05, though this number was arbi- 
trary and aside from an obvious rescaling, does not affect the 
results obtained. 

To observe behavior in the limit of metastability as dis- 
cussed in the previous section, the rate of unbinding must be 
set to be small compared to the inverse metastable equilibra- 
tion time in the bound state. Hence we chose the unbinding 
rate co = 10~ 5 and c\ = 5 x 10~ 4 . 

As the simulation was running, the tether point was moved 
at velocity v which was typically small, The average spring 
force (/) was measured as a function of v for a given spring 
constant k. The results of a run of 3 x 10 10 steps are shown in 
Fig. Eta). 
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FIG. 8. (a) The force and power versus velocity for the one dimen- 
sional model described in the text. Here the spring coefficient k — 15. 
The line going through the power is a cubic spline fit used to more 
accurately determine the maximum value of the power. The velocity 
and power are both in units of 10~ 6 . (b) A plot of the efficiency of 
the motor as a function of the spring constant k. Three separate runs, 
each of 3 x 10 10 steps were used to determine the error bars, for each 
data point shown. 



We are interested in the limit of large k although this is hard 
to achieve numerically owing to exponentially long relaxation 
times. The parameters used allow us to probe up to k = 20. In 
the large k limit, the minimum value of V„„-„ needed to bind, 
is kl 2 /2, where I is the maximum amount the spring will need 
to stretch from the tether to the binding site. Because 8 = 0.1 
and L = 1, this implies / = 0.9. Therefore our formula for 
the efficiency, given this input energy is tj = P max /{ckl 2 /2), 
where P max is the maximum of the power versus velocity, as 
shown in Fig. [Hfa). Plotting the efficiency for different values 
of k ranging from 5 to 20 yields the points in Fig. 13b). De- 
spite the fact that the relaxation time for metastable relaxation 
of the system varies over more than two orders of magnitude, 
the overall efficiency is almost constant. We expect at higher 
value of k, the efficiency will eventually drop owing to the 



fact that the average unbinding time becoming smaller than 
the metastable relaxation time. 

What the above analysis shows is that in the limit where 
the photon cross section can be made arbitrarily low, the ef- 
ficiency can be adjusted to be constant, independent of the 
photon energy. This is accomplished by choosing a large 
spring coefficient. In reality with photon energy of 2eV and 
kgT ss 1 /40eV, photon flux would have to be far too low for 
this optimal regime to be realizable. However the above anal- 
ysis also shows that the efficiency can be substantially in- 
creased by choosing larger k at the expense of lowering the 
cross section. The energy delivered in one cycle should scale 
as E c = kL 2 . This can be increased to be substantially larger 
than thermal energies but at the expense of a long relaxation 
time proportional to exp(E c /ksT). As noted earlier, it should 
be possible to make E c about IksT, with reasonable parameter 
estimates. 



VIII. CONCLUSIONS 

Here we have analyzed the viability of converting photons 
to mechanical energy using a device composed of an canted 
polymer brush tethered to a lower plate but able to bind its 
other ends to sites on an upper plate. Photons can dissociate 
these ends from binding sites. By a combination of analytical 
and numerical arguments we showed that in steady state, this 
produces net mechanical power. 

The system is inspired by biological motors such as myosin 
II that bind to actin and is dissociated by the binding of ATP. 
The analysis used here could also be applied to such systems, 
however in reality they contain many more stages. In general, 
these kind of systems are classified as "thermal ratchets" 19 , 
where the system can be thought of as moving in a wash- 
board potential in the presence of thermal noise. Though that 
description can be very useful in understanding the general 
principles behind the operation of such motors, in the present 
case, we are trying to model the system in more detail than 
such models can afford. Instead we have described the sys- 
tem using Langevin dynamics and also two coupled Fokker- 
Planck equations similar in spirit but not identical to previous 
approaches 1415 . The difference here is that for potentials to be 
a sensible model for a polymer tethered to a single point, they 
cannot be periodic. Such modelling allows us to see how vary- 
ing microscopic parameters affect the power and force char- 
acteristics. 

The analytical results on the Unbound Equilibration Model 
and extensive one dimensional simulations, show that the 
force applied by the device can be made arbitrarily large at 
the expense of having exponentially long relaxation times. At 
a given photon flux, the production of large forces from sin- 
gle polymers imply the need for a very low cross section of 
interaction between the photon and the bound end plus bind- 
ing site, as long relaxation times are required. Therefore there 
is a trade off between this force and the speed the device can 
move. This may be circumvented to some extent by stacking 
transparent devices of this kind, so that even though the cross- 
section of interaction of an individual photon with a given 



10 



layer is low, it will eventually be captured by one layer. 

Therefore it appears that there is no theoretical obstacle 
to prevent the photo-mechanical conversion of energy in this 
manner, however it represents a significant experimental chal- 
lenge. 

It has been pointed out 22 that there is an important distinc- 
tion between artificial molecular switches and artificial molec- 
ular machines, the former being a fraction of a penny, and the 
latter being extremely challenging to create. The approach 
investigated here is closer to biological motors than other pro- 
posals, and should be more forgiving about randomness, ei- 
ther during fabrication, or due to thermal motion, than ap- 
proaches that require precise chemical synthesis of molecules 
capable of sophisticated conformational changes 6 . However 
its experimental realization is still quite clearly a formidable 
task. 
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